Astronomy & Astrophysics manuscript no. paper' AA-201 1-17151 '2arxiv 
February 27, 2013 



©ESO2013 



Rotational splittings for slow to moderate rotators: 

Latitudinal dependency or higher order effects in O? 

R-M. Ouazzani 1 , and M-J. Goupil 1 



LESIA, UMR8109, Universite Pierre et Marie Curie, Universite Denis Diderot, Observatoire de Paris, 92195 Meudon Cedex, France 
e-mail: rhita-maria.ouazzani@obspm.fr 



o 



O 
CN 

P? 
CO 

6 



> 

oo 
m 
^r 

en 

o 

CN 



X 



Received April 1, 2011; accepted May 16, 2011 



ABSTRACT 



Context. The unprecedented photometric quality reached by the CoRoT and Kepler space missions opens new prospects for studying 
stellar rotation. Information about the rotation rate is contained on the one hand in the low frequency part of power spectra, where 
signatures of nonuniform surface rotation are expected, and on the second hand in the frequency splittings induced by the internal 
rotation rate. 

Aims. We wish to figure out whether the differences between the seismic rotation period as determined by a mean rotational splitting, 
and the rotation period measured from the low frequency peak in the Fourier spectrum - observed for some of CoRoT's targets - can 
provide constraints on the rotation profile. 

Methods. For uniform moderate rotators,perturbative corrections to second order and third order in terms of the rotation angular 
velocity Q., must not be neglected. These effects, in particular, may mimic differential rotation. We apply our perturbation method 
to evaluate mode frequencies accurate up to Q. 3 for uniform rotation. Effects of latitudinal dependence are calculated in the linear 
approximation. Numerical results were obtained for selected models of the upper and lower main sequence. For the latitudinal depen- 
dence, we adopted two types of rotation profile: one with rotation uniform in depth, and one with a solar-like tachocline. 
Results. In models of f3 Cephei pulsators, upper main sequence stars, third order effects become comparable to that of a horizontal 
shear similar to the solar one at rotation rates well below the breakup values. We show how a clean signature of the latitudinal shear 
may be extracted. Our models of two CoRoT target HD 181906 and HD 181420, which are solar-like pulsators, represent lower main 
sequence objects. These are slow rotators and nonlinear effects in splittings are accordingly small. We use data on one low frequency 
peak and one splitting of a dipolar mode to constrain the rotation profile in HD 181420 and HD 181906. 

Conclusions. The relative influence of the two effects strongly depends on the type of the oscillation modes at stake in the star and on 
the magnitude of the rotation rate. Given mean rotational splitting and the frequency of a spot signature, it is possible to distinguish 
between the two hypothesis, and in the case of differential rotation in latitude, we propose a method to determine the type of rotation 
profile and a range of values for the shear. 

Key words, stellar oscillation - stellar rotation - seismology - differential rotation 



1. Introduction 

The CoRoT and Kepler space-borne missions with their unin- 
terrupted observations spanning a long-time interval promise a 
wealth of data suitable for studying stellar rotation. Information 
about internal rotation is contained in characteristic spacings 
(known as the rotational splittings) which appear in the power 
spectra of light curves. Mean values of the rotational splitting 
have already been de termined for a number of stars -HD 181420 
dBarban et all 120091) . HP 49933 dBenomar et a ll 120091 V1449 
Aql dBelkacem et all 120091) and HD 181906 JGarciaet all 
l2009h - observed by CoRoT. These values yield some informa- 
tion about average rotation rates sampled by modes detected 
in these objects. Since all the data concern p-modes, the mean 
values mostly reflect the rotation rate in the outer layers. To 
probe deeper layers, we need splittings for gravity modes. In 
three Fourier spectra analyzed so far (HD 181420, HD 181906 
and VI 449 Aql), low frequency peaks were found but were at- 
tributed to the effects of spots on the rotati ng stellar surfaces 
(for a complet e review on spot modeling, see ICollier CameronL 
2002; M osser et all l2009h . Rotation periods deduced this way 
were found to be different from those determined from the split- 
tings. This is not surprising. The spots only give access to the 



surface rotation rate at the latitude of their location, whereas the 
splittings yield a mode-dependent mean value of the interior pro- 
file. 

The linear relation between the splittings and the rotation 
rate, £2, follows from the first order p erturbative treatment of the 
Coriolis acceleration (lLedouxlll945l) . At moderate rotation rates, 
the perturbative formalism may still be applicable but we have 
to go beyond the fir s t order (iReese et all 120061 rOuazzani et all 



2009t ISuarez et all 120101: iBurke et all 1201 ll). ISaiol (119811) 



Gough & Thompson! dl9901) and lDziembowski & Goodd dl992h 



(hereafter DG92) derived oscillation frequencies including sec- 
ond order corrections in Q. Such corrections arise from the 
higher order effects of the Coriolis acceleration and the lowest 
order effects of the centrifugal acceleration, which causes distor- 
tion of the stellar structure. The oscillation frequencies no longer 
depend on CI in a linear way but linear rotational splittings may 
still be recovered from frequency differences between prograde 
and retrograde modes of the same degree and order. This sim- 
ple property is lost when rotation couples modes with close fre- 
quencies. The formalism must then be modified (DG92, Soufi 
et al. 1998) and the recovery is more difficult. The cubic effects 
in Q make it still more complicated. The resulting difference be- 
tween prograde and retrograde modes becomes dependent on the 
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mode's azimuthal order, m, in such a way that it may be misin- 
terpreted as the effects of a latitudinal dependence of the rotation 
rate. 

The question arises whether in the presence of significant 
nonlinear effects it is still possible to extract the values of lin- 
ear splittings, which provide integral constraints on differential 
rotation in the interior. The next question which we ask in this 
paper is what may be learnt by combining such constraints with 
data on low frequency peaks that are attributed to spots and yield 
information on the surface rotation rate. We expect that the an- 
swers depend on the type of pulsator, and the characteristics of 
its observed modes. Here, we specifically consider two very dif- 
ferent types of main sequence objects, /3 Cephei and a solar-like 
pulsators. In the first case, we are dealing with a massive star 
where low order p- and g-modes are unstable, and in the second 
case with low-mass stars where high order p-modes are damped 
but stochastically driven. 

The paper is organized as follows. Sect. 2 gives the basic 
theoretical framework for this study. In Sect. 3, we focus on the 
effects of cubic order and near degeneracy contributions to pul- 
sation frequencies and present numerical results for a selected 
model of a /? Cephei star. Explicit expressions for the rotational 
splitting in the case of latitude dependent rotation profile are 
given in Sect. 4. Departures from a linear dependence on m are 
then compared with those connected to the nonlinear effects in 
Q. Prospects for disentangling these two effects are discussed in 
Sect. 5, using the case of a/3 Cephei star as an example. Sect. 6 
is devoted to two CoRoT targets, the solar type stars HD 181906 
and HD 18 1420, for which we combine the splitting with the low 
frequency peak and make some inference on the rotation pro- 
files. Sect. 7 is dedicated to conclusions. 



2. Perturbational treatment of uniform rotation: 
effects on pulsation frequencies 

In the presence of rotation, the centrifugal and Coriolis acceler- 
ations come into play. The centrifugal force affects the structure 
of the star and distorts its shape. The resonant cavity is changed 
and with that the oscillation frequencies. The Coriolis force en- 
ters the equation of motion and affects the motion of the waves, 
hence the frequencies of normal modes. As rotation breaks the 
spherical symmetry, it lifts the frequency degeneracy, introduc- 
ing a dependence on the azimuthal order, m. Without rotation, 
mode frequencies only depend on the radial order, «, and the 
angular degree, (, and are 2£ + 1-fold degenerated. 

To first order in the rotation rate, SI, the normal mode fre- 
quencies in the inertial frame are given by co„ j Jn = a» n ,f,o + 
m Q/3,,/ . The explicit form of the last coefficient is provided in 
Eq.( IB.5b . Here, w e only note tha t C,,/ = 1 - fl n ,c is known as the 
Ledoux constant dLedouxlll945h . 

Characteristic spacings appear in the spectrum, such as the 
rotational splitting that we define here as: 



%i,t,m - W)i,{,m - Wn,t,-m 1 2m 



(1) 



This gives us a basis for determining the rotation rate. From now 
on, we drop the subscripts (n, (, m) for to and the splitting, unless 
there is an ambiguity. We also use <x = co/Slk (where Slk = 
yGM/P? is the break-up frequency) instead of u>. 



2.1. Equilibrium configuration 

We here consider the case of uniform rotation. The stationary 
equation of motion in an inertial frame of reference is: 



VP 

(v -V)vo = V0 = F 

P 



(2) 



Where in the left hand side, vo = SI Ar = SIqi sin^e^ in the 
spherical basis is the velocity field due to rot ation at the angula r 
velocity SIq, 6 being the colatitude (see e.g. lUnno et all I1989I) . 
P, p and <p are the pressure, density and gravitational potential, 
respectively. For a rotating star, the left hand side corresponds to 
the centrifugal acceleration F = -SI x (SI x r), whose effect on 
the equilibrium structure is twofold : on one hand, a spherically 
symmetric perturbation, which mainly modifies the gravity; on 
the other hand, ^-dependent perturbations, which is responsible 
for oblateness. Then all the equilibrium quantities, X, are well 
approximated by: 



X(r,0) - X(r) + X 22 (r)P 2 (cos#), 



(3) 



The spherically symmet ric part is t hen obtained by (see for ex- 
ample Ki ppenhahn & WeigertL Il994h 

dp . GM r 2 2 

T = -PEeff. where g eff = — -rfl 

dr r z 3 

The non-spherically symmetric part is obtained by (see DG92) 

(4) 
\ r- J/ 

prSl 2 (d\np\((f>22 „ 2 1 



P22 = -pr SI \ — Sl + -I, 



P22 = 



g \dlnr/\r 2 3 1 



where 



1 d / 2 d<fe\ 6 

?*** - r 2fe=4^ Gp22 . 



(5) 



(6) 



The boundary conditions can be found in ISoufi et all (1 19981) 
(hereafter S98). 

2.2. Oscillation frequencies up to cubic order in SI 

Using expansions of the type given in Eq.Q for the oscillation 
quantities, i.e., p' = p' + p' 2 , the oscillation system is then ex- 
pande d up to the cubic order. According to S98 (see also lKaramiL 
J2Q08!), the oscillation equation then becomes: 



u 



(A + eB)£ + e 2 (D + eC)£ + 0(e A ) = 



(7) 



e being equal to Sl/Slt, where Slk - -y/GM/R 3 is the break-up 
frequency. The operator A represents the basic linear oscillation 
operator including the spherically symmetric perturbation due 
to rotation (through g e fr). The operators B and D, respectively, 
contain the effects of the Coriolis force and of non-spherically 
symmetric distortion. The operator C shows that a coupling be- 
tween the non-spherically symmetric distortion and the Coriolis 
force exists. 

Like in S98 and in lKaramil J2008I) . parts of the Coriolis and cen- 
trifugal distortion effects are included into the pseudo-zero th or- 
der eigenvalue system. This way, we are able to solve the eigen- 
value problem up to cubic order without having to solve the 
successive equations for the eigenfunctions at each order. The 
solution yields eigenfrequencies ctq, which include parts of the 
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frequency shifts induced by rotation. To single out various con- 
tributions and emphasize the m-dependence, we write it in the 
form 



10) 



.mi 



O"0,m = C" + 0~ 2 + <T\ 



eisen eieen 

+ (T^ + CT, B 

2,/n 5.m 



(8) 



CT* = cr,„ ± 



V A ». + ^ 



(13) 



where cr,„ = - and A,„ = 



where cr^ is the classical zero* order frequency ignoring all 

effects of rotation, <t 2 ' is the correction resulting from the 
spherically-symmetric part of the centrifugal distortion, and the 
next three terms give the contributions of consecutive orders 
in Q resulting from the Coriolis acceleration. The linear term, 
cr lm = milfXlk, is complete. The higher order terms only in- 
clude the parts resulting from the poloidal component of £. The 
remaining contributions to frequency shifts up to 0(Q. 3 ) are cal- 
culated as integrals involving the eigenvectors £q ( see S98). To 
this accuracy, the complete expression for eigenfrequencies in 
uniformly rotating star is given by 



0"m — @~0,m ' @~c,m 

with o- cjn 



T D T D C 

<rL + f*L + <*%„ + of + o%_, 



(9) 
(10) 



where the exponent T marks contributions from the Coriolis 
force acting on the toroidal component of £, the exponent D 
those arising from the non-spherically symmetric distortion, and 
c those resulting from coupling of the two effects. 



where the sign +(resp.-) corresponds to mode k\ (resp. £3). Note 
that if \ < H,„\ « |A,„|, the effects of coupling are small and mode 
frequencies are approximately described by Eq.(|9}- 

3. Relative magnitude of the different contributions 
for uniform rotators 

In this section, we adopt as an example a stellar model with a 
simple equilibrium structure for which uniform rotation in depth 
is assumed. All the results presented in this Section have been 
performed for an 8 .5M ZAMS stellar model described in Table 
Q] rotating at Q - 15% Q* (v n = RQ = gSkm.s 1 ) which is 
representative for this type of star (IStankov & Handler! l2005h . 

Table 1. Stellar parameters of the ZAMS model (Sect.2 and 3.) 



M = 8.5 M 


R = 3.96 R Q 


L = 13 x 10 3 L o 


X = 0.7 


P c = 8.8 x 10 3 dyn.cnT 2 


p c = 3.6 x ltr 9 g.cnr 3 


o- rot = i5%n A 


v rot = 95km. s _1 



2.3. Near degeneracy 

The standard perturbation approach is invalid if rotation couples 
modes with close frequencies. Treatment of such cases requires 
modification, which in the context of stellar pulsations was first 
used by Chandrasekhar & Lebovitz (1962) and developed later 
by DG92. In the case of latitude-independent rotation profiles, 
only modes with the same m and £s of the same parity are cou- 
pled. In the present work, we study the coupling of two reso- 
nant modes denoted kj for (rij, (j, m), j — 1,2 with frequencies 
<T\ > o~2- For the range of stellar models we are interested in, 
calculations reveal that near degeneracy occurs for quite a large 
number of modes. 

Near degeneracy is taken into account by searching for solu- 
tions of Eq.(|7]i in the form 



£=£^o,*, +& 7=1,2 



Stellar models and adiabatic oscillation frequencies are 
computed with the evol ution code CESAM2k dMorell Il997t 
iMorel & Lebretonl 120081) and the WarM (Warsaw Meudon) os- 
cillation code, respectively (see S98). All numerical results pre- 
sented here and in the rest of this paper concern dipolar (( = 1) 
modes. 



3.1. Frequency contributions up to cubic order 



and cr,^ 

3,m 



Here we quantify the implicit contributions (o" 2 
from Eq.®), and compare them with the corrective terms of the 
same orders given in Eq.dTOb. To extract these implicit contri- 
butions from the pseudo-zero* order eigenfrequencies we make 
use of their symmetry properties in m. Thus, we get 



0% 



(O"0,m + C"o - m ) 



-CT, 



(0, 



(14) 



with 



(11) and 



htkj 



k= l,N, 



where the eigenfunction correction £ c is composed of all non- 
resonant modes. The standard procedure leads to a linear system 
of equations for 9\i and the following condition for a non-zero 
solution: 



(cr k] - crt s )(o- kl - (Tm S ) 



w- 







(12) 



where the frequencies cr Al and cr>, are given by Eq.©. The cou- 
pling term 7Y„, corresponds to integrals containing second and 
third order contributions (see S98 for more details). The solu- 
tions of Eq.dTSt. denoted by cr deg , provide the desired eigenfre- 
quencies: 



3,m 



(f"0,m - 0"O,-m) £2 
- m p 

2 af 



(15) 



Figure Q] shows all the contributions of different orders to the 
rotational frequency corrections for g-modes (gi4 to gi) and p- 
modes (pi to pig). We note the significantly different pattern for 
p- and g-modes. In the latter case, the second order correction 
due to the Coriolis acceleration dominates the rest. It grows lin- 
early wi th the radial order (in absolute value) as was already 
stated in lBallot et alj d2010l) for polytropic models. The Coriolis 
effect remains dominant at the third order, whereas distortion 
seems to have no effect. The coupling between the Coriolis ac- 
celaration and distortion (o"3 C ), is negative, thereby reducing the 
impact of the former. 

For p-modes, the dominant second order term is cn£, as could 
be expected for modes which are mostly confined in superficial 
layers where the the role of centrifugal acceleration is highest. 
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Fig. 1. Contributions of different approximation orders to / = 1 , m — mode frequencies: <X2 e igen (o"3eigen) represent the implicit 2 nd 
(3 rd respectively) order contribution to the eigenfrequency; ct 2 t and <x 3T denote 2 nd and 3 rd order frequency corrections due to the 
Coriolis acceleration; ct 2 d and <x 3D represent 2 nd and 3 rd order frequency corrections due to the centrifugal distorsion; cr^c comes 
from the coupling of distorsion and Coriolis effects. The contributions are plotted as a function of the radial order n for an 8.5 M Q 
ZAMS model rotating uniformly at 15% Qa, i.e. around 95 km s~ : (see stellar parameters in TablefT}. The frequencies are scaled by 
£lk - \/GM/R 3 : <x = co/Qk. We recall the use of a negative radial order n for gravity modes and a positive one for acoustic modes. 



In this case too, the dominan t term gro ws linearly with the ra- 
dial order, as noted in lGoupill(l2009l) and lReese et al.1 (l2006h for 
polytropic models. At the third order, <x" sen and <x° are of the 
same magnitude but with opposite signs and cancel each other 
to some extent. 

To sum up, as expected, for g-modes the most important con- 
tributions are related to the effects of the Coriolis acceleration, 
whereas for p-modes, we must take into account both the im- 
plicit eigenfrequency terms (due to the part of Coriolis force in- 
cluded in the pseudo-zero* eigen-system) and the effects of cen- 
trifugal distortion. In Table lD.3l numerical values of the different 
contributions are listed. 

3.2. Near degeneracy corrections 

Here we consider the same sequence of € = 1 modes as in the 
previous section but we now take into account the coupling of 
each mode with the nearest t = 3 partner. The coupled pairs 
must be of the same azimuthal order m. We use Eq.(fT3l here to 
calculate the frequency shift caused by such a coupling. 

In order to compare the magnitude of this near degeneracy 
effect with the second and third order contributions shown in 
FigQ] we depict in Fig[2]the frequency differences between com- 
putations with and without near degeneracy being accounted for. 
Figure [2] shows that near degeneracy primarily affects the p- 
modes (n > 0). This is not surprising as they are more sensi- 
tive to the outer regions which are more affected by distortion, 
a dominant factor in the coupling coefficient ( H. Moreover, this 
correction is found to be of the same magnitude as the other 
second order corrections (see FigJT} but of the opposite sign. 



1 = 1 mode, m = 1 i+ 

m = —-X — - 

m = - 1 — -&K- 









...... 



.-./ 



"A A ■■ VN. A - A' A 1 A A PSTTTf ■ r ■ v 
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3 - 

2 - 

1 - 

-3IOK- 

-1 ■ ' 

-16-14-12-10-8 -6-4-2 2 4 6 8 10 12 14 16 
radial order n 



Fig. 2. Frequency differences between degenerate c^^ and non- 
degenerate cr nodeg solutions scaled by fi A = V GM /R 3 (°" = 
tL>/£lk) as a function of the radial order, n, for 6 - 1, m - -1, 0, 1 
modes. The computations were made for the same model as in 
Fig. O. 



Hence the overall effect of the distortion is reduced. However, 
this is not a universal property. As we may see in Eq. rtOV the 
coupling always causes an increase of frequency separation be- 
tween modes but the sign of the shift is mode-dependent. In any 
case, rotational mode coupling is an important effect, especially 
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for p- modes. Taking it into account (as shown in ISuarez et all 
120101) extends the validity domain of perturbative methods. 

Finally, we notice that the sectorial components of the I — 1 
triplet are modified by roughly the same amount which implies 
that the rotational splitting should not be strongly affected by 
near degeneracy (see Sect. 14. 2t . This is expected because in this 
case it enters as a third order effect. 

4. Rotational splitting for uniform rotators 

The rotational splitting can be defined as : S m = (cr m - cro)/m. 
One also uses S m = <x m - <x m _i. In this work, we use a scaled 
expression of the rotational splitting (Eq. ([!}): 



S m = cr m -cr_ m /2m 



(16) 



These various definitions are equivalent only at first order in the 
rotation rate, Q, and equal to the linear splitting: 

S «>= 7V f f K m (T,e)Cl(T,0)d66r (17) 

"A Jo Jo 

where the ana lytical expression for the kernels K m is given in 
iGoupill (1201 ll) and references therein. At higher orders in terms 
of £1, the two first definitions are contaminated by the effect of 
asphericity, which introduces an antisymmetric component in 
the frequency as a function of m. We choose to remove this sec- 
ond order contribution using the splitting expressed in Eq.dToTi. 



4.1. Cubic order effects on the splitting 

Assuming a uniform rotation Q = Qo, the splitting including 
frequency correction due to cubic order effects is given by: 



i cubic 



LIq ±1q I Qi 



L 3,m 



Ll k cr \ Ll k J 

+ TL, + T?„ + T 



3,m 



L 3,m 



3,m 



L 3,m 



(18) 
(19) 



where <x = (Oo/Qk 1S the normalized frequency of the corre- 
sponding axisymmetric mode, Qo is the uniform rotation rate, 
and j3 is the integral of the kernel K m over the star (see Appendix 
IB Eq.dUl). 

T3 m contains the implicit third order contribution as well as 
effects due to the Coriolis acceleration, the distorsion, and the 
coupling of the two. From now on, we define the departure from 
a linear splitting as follows: 



c C cubic _ C cubic / Q o\o ^o/^o\ T 



4.2. Near degeneracy correction 



(20) 



According to the formalism explained in Section 123] if we con- 
sider the coupling of t — 1 and I — 3 modes, let the degenerate 
frequency of ( — 1 modes be: 



t=\,i 



07=1. 



07=i 



p-^ + ^ + <Hl (21) 



m L -^m 



A,„ + 



where cr f=1 m and crc=xm are non degenerate frequencies given by 
Eq.©, and A,„ is defined in Eq.(fT3l. Then the splitting account- 
ing for near degeneracy is given by: 



Q deg _qND 1 /a 



A_i) (22) 



where 



tND 



(T\ +1 - <T\ -1 



(23) 



with ND standing for non-degenerate. Note that S t f D . , con- 

& & t=l,m=±l 

tains cubic order contributions mentionned in the previous sec- 
tion. The contribution of near degeneracy to the splitting is then 
givenbyS^-S-^. 

It is worth recalling that neglecting all cubic order contri- 
butions in Eqs. (O, ( fTUb and ( fT3b results in Ai = A_i and 
fii = < H-\. In that particular case, degeneracy does not con- 
tribute to the rotational splitting, and the rotational splitting is 
linear in Q up to second order and satisfies : 



S, 






(24) 



4.3. Sensitivity to the nature of the eigenmode 

Figure [3] displays the near degenerate contributions to the split- 
ting for p- and g-modes together with the cubic ones for a ZAMS 
stellar model (Table [TJ. The values of these different contribu- 
tions are given mode by mode in Table lD.2l 

In Fig[3] (left), the Coriolis correction T^ dominates for g- 
modes (by roughly a factor 10 2 over other contributions) and 
decreases with the radial order to a roughly constant value for 
low \n\. The scale is too large in this figure to see the behavior of 
T" se ", T^, T^ and the near degeneracy contribution to the split- 
ting, but we refer to Table ID. fl where it is shown that there is no 
asymptotic behavior for these four contributions. Near degener- 
acy is fully negligible for all g-modes except for the n — -1 
mode that actually is a mixed mode. The g-mode spectrum is 
much denser than the p-mode one but as shown in Fig. IC.2l in 
Appendix ICl the coupling term H m is much smaller than A m . 
This is due to the fact that distorsion effects are small for g- 
modes and therefore the overall (second and third order) contri- 
bution to H„, remains small. 

In order to emphasize possible asymptotic behavior, the con- 
tributions to p-modes have been divided by cr 2 ,. In Fig. [3] (right), 
the implicit cubic order and the near-degenerate contributions 
dominate for p-modes but with opposite signs and therefore 
roughly compensate each other. Hence, as was the case for pul- 
sation frequencies, the near degeneracy correction tends to re- 
duce over-estimated contributions to rotational splittings. The 
contributions dominated by centrifugal distorsion, and in par- 
ticular the near degeneracy one, scale as cr 2 . In AppendixlAl it is 
shown that near degeneracy contributes to the rotational splitting 
only if third order effects are taken into account (see Eq. dA.5l i). 
However it is also shown that second order effects CH{) - dom- 
inated by distorsion CH2) for P modes - are involved. This ex- 
plains why the near degeneracy frequency variation follows a <x 2 
behavior for p-modes. 

Similar conclusions for more evolved models with more 
complex structures are found for pure p-modes and pure g- 
modes. However these complex structures also give rise to mixed 
modes for which all effects contribute equally, and a precise in- 
vestigation, mode by mode, has to be done for each equilibrium 
model. This will be investigated in details in Sect. [6] 
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5. Effects of latitudinal shear on the splitting 



lHansen et al.l d 1977b derived the expression for the rotational 
splitting of adiabatic nonradial oscillations for slow differ- 
ential (steady, axially symmetric) rotation Q.(r,9) and ap- 
plied it to numerical models of white dwarfs and of massive 
main sequence stars assuming a cylindrically symmetric ro- 
tation law. In the solar case, the effects of latitudinal differ- 
ential rotation on theoretic a l frequencies were investigate d by 
Gough & Thompson! (1 1990l).lDziembowski & Goodel(ll99ll) and 



Dziembowski & Gooddd 19921) . 
In order to be able to compute the splittings from Eq.dTTt. one 
must specify a rotation law. It is convenient to assume the fol- 
lowing form: 



-'max 

Q(r, 8) = ^ ^2s(r) (cos 9) 2 



(25) 



where 6 is the colatitude. The surface rotation rate at the equator 
is Q(r = R, 9 = n/2) = Q (r = R). 

Note that in the solar case, Q2 and Q4 are negative and the equa- 
tor rotates faster than the poles. 



Inserting Eq.(f25b into Eq.(fT7b yields the following expression 
for the splitting: 



1 ^R 1 s=2 

S m = — I flo(r) K(r) dr + — V m 2s H m , s (Q) (26) 

"* Jo "k ^Tq 



The expression for H ms (Q) can be found in AppendixlBl 



5. 7 . Latitudinally differential rotation only 0,(6) 

In this case, for which r otation is assum ed to be uniform in depth, 
the splitting becomes: (iGoupilLbOllh 



Q 1 

Q * Q * to 
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Expressions for R s and Q s can be found in AppendixlBl /3 and y 
depend on the radial and horizontal components of the mode: 



= -1/1 £, 

Jo 



y= -i/i I ^pordr 
Jo 

I being the inertia of the mode: 

I 



f te 2 +A4 2 ) P0 r 2 dr 
Jo 



(28) 



(29) 



(30) 



For the sake of simplicity, we restrict our study to s = 1 in 
Eq.GoY The rotation law can then be expressed as: 

Q(0) = Q o -AQcos 2 # (31) 

Where Q = Q(r = R, 9 = n/2), Q 2 = -AQ and Q 4 = 0. 
After some calculations provided in Appendix IB] we are able to 
express the splittings of several ({, m) combinations: 

a. (€, m) = (1, ±1) modes: 



b. (^m) = (2, 



Mat 



clat 

:1) modes: 

Qo 



=2,m=l 



Q,t 5 Qo 



1 AQ 
S+ --(8 r -3^) 



(32) 



(33) 



(Am) = (2, ±2) modes: 

_ Qo 

' Q* 

In the solar case, /3 ~ 1 and /? >> \y\ for the excited high fre 
quency p-modes. Then the (I - 1, m = 1) splitting is: 



Mat 
V=2,i 



1 AQ 

2/?---(2 r+/3 ) 



(34) 



clat ^ — " 



1 



1 + - (Q 2 /Q 



3/7 Q4/Q0) 



(35) 



Qo 

Qa 

With Q 2 /Qo = -0.127, Q4/Q0 = -0.159 , one obtains a depar- 
ture of |Si/Qo-l| - 0.04 from a linear splitting, i.e. a 4% change 
in the solar case. 

For upper main sequence stars, excited modes are around 
the fundamental radial mode and may be mixed modes with 
|/8| ~ \y\ ~ 1/2. This leads for instance to |S U /Q - 1| ~ 1% for 
Q2/Q0 and Q4/Q0 equal to a half of the solar values. We recall 
that for the stars treated in this article, we have taken Q2 = -AQ 
and Q4 = 0. 

As shown in Eqs. ( TZTT l to Eq.(t30b. the additional term due 
to the latitudinal shear strongly depends on the eigenfunction of 
the mode, through the radial and horizontal components of the 
displacement. Therefore, we investigate the scaled contributions 
for different types of modes. Again, the plots presented in Fig|4] 
have been performed for an 8.5 M , 4 R Q ZAMS model, rotating 
at Q si 15% Qa (vp_ = RQ = 95km.s~'). Once again, com- 
putations for more evolved stellar models give similar results - 
except for mixed modes around n = 1 - even with stronger rota- 
tion rates. 

Figure [4] shows the dependency of the /3 and y integrals on 
the p or g nature of the mode. Globally, for g-modes, at a given 
rotation rate, we expect a small contribution of the latitudinal 
shear to the splittings, smaller for 5i ( i than for S2J and S2,2, 
whereas for p-modes, the contribution is quite important for S24. 

For the sake of simplicity we will focus on £ — 1 splittings 
from now on. Note that the investigations presented in the 2 next 
Sections have been addressed for 1-2 splittings as well and had 
lead to the same conclusions. 



5.2. A tachocline-like rotation profile: Q(r, 9) 

Let us refine the modeling of the rotation profile, assuming a 
rotation profile in depth as in the solar case. Rotation is uniform 
(Q = Qo) in the inner layers and differential latitudinally- as 
expressed in d3Tt - in the convective envelope. r cz being the radius 
of the inner stable layer: 

for r < r cz , Q(r, 9) = Q 

r > r cz , Q(r, 6) = Q - AQ cos 2 9 (36) 

Therefore, after some calculations similar to those presented in 
Appendix iBl equation (HTt is no longer relevant, and should be 
replaced by: 



Iat _ ^0 / AQ 1 

S( = l - m = l - OT r " Qo" 5^ : 



(37) 



1 f R 
where p cz = y J [g - 2£& + (A - l)4 2 ]p r 2 dr (38) 

Note that for p modes for instance, /3 CZ is smaller than /3 in 
the whole star - /3 CZ ~ 0.45 whereas {3 ~ 1 -. Therefore in the 
case of a tachocline-like profile, the effect of a same latitudinal 
shear AQ is smaller than in the case where rotation is uniform in 
depth. 



6. The case of a/3 Cephei on the main sequence 

Massive stars on the main sequence are usually fast rotators and 
their fast rotation affects their internal structure as well as their 
evolution. Rotation of /3 Cephei stars extends from slow rota- 
tional velocity (v < 50 km/s) to extremely rapid ones (v > 250 
km/s) (IStankov & Handler! 120051) . These stars usually have a ra- 
diative envelope which may or may not be in latitudinal differ- 
ential rotation. For these fast rotators, one can wonder whether 
cubic order or near degeneracy contributions dominate over the 
effects of latitudinal shear. 

Here we investigate the importance of deviations from a lin- 
ear splitting for an evolved main sequence model of an 8.5 M 
star with a 5.07 R radius (see Table|2]for the stellar parameters 
of the model). 



Table 2. Stellar parameters of the /3 Cephei model in Sect. 5 



M = 8.50 M G 


R = 5.07 R G 


L = 2 x 10 37 erg.s^ 1 


age = 20 My 


X = 0.713 


Z = 0.014 


a = 1.76 





6. 1 . Departure from a linear splitting as a function of the 
frequency 

According to the region to which a mode belongs, very different 
types of behaviour are observed for the different contributions to 
the splitting, as illustrated in Figj5] 

Figure [5] displays the departure from a linear splitting due 
to cubic order effects with and without near degeneracy correc- 
tions, and to latitudinal differential rotation (uniform in depth). 
It has been computed for £ - I modes with radial orders ranging 
from -10 to 5. The parameters of the stellar model considered 
here are given Table [2] It is rotating at 20%Qt (left) and 30%Qt 
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(right). On these plots we observe three regions where the be- 
havior of 6S c,m differs. 

In the g-mode region (« < -2), as the cubic order effect is 
proportional to (^-J Jr (see Eq.d20"li), the smaller the frequency 
is, the higher cubic order terms are, and this effect increases with 
the rotation rate. In this region, cubic order effects overtake those 
from latitudinally differential rotation. 

In the p-mode region (n > 2), latitudinal shear effects are 
larger than for g-modes. Cubic order effects are of same order 
of magnitude as the contribution from latitudinally differential 
rotation with a shear of only 12.7 %. 

In between, in the mixed mode region, cubic order effects 
with degeneracy corrections are of same order of magnitude as 
effects from latitudinally differential rotation for all the consid- 
ered shears. 

In order to get a better understanding of these different types 
of behavior, one can look at Fig lC.ll in Appendix [C] Figure 
IC. II displays mode inertia (defined in Eq.(l30"il) as a function of 
the radial order, for the model described in Table |2] uniformly 
rotating at 30% fit . From this figure, the three frequency do- 
mains related to the nature of modes are clearly visible: below 
g3 (n e [-10, -3]) are pure g-modes, above p2 are pure p-modes, 
and in between is located the mixed mode region. 



Figure IC.2I which has been computed for the model de- 
scribed in Table [2] uniformly rotating at 30 % O^, shows that 
near degeneracy has a very small effect on g-mode splittings: 
although frequency differences (A ; „) between I - \ and I — 3 
modes are small for g-modes (the y axis is in a logarithmic 
scale), the fact that "H,,, (which quantifies near degeneracy) is 
dominated by second order terms (due to distorsion) and takes 
small values for g-modes makes near degeneracy more or less 
negligible for g-modes. In the mixed mode region, although "K,„ 
increases, the frequency differences between I - \ and € = 3 
modes A,„ also increase. Consequently, near degeneracy should 
not be important in the range either. Once again, as mixed modes 
are very sensitive to the evolutionary stage, it is diffucult to make 
general statement for all the stellar models. 



In the p-mode domain, A,„ decreases, while "7Y„, increases 
causing the near degeneracy contribution to the splitting to in- 
crease. Therefore, the near degeneracy effect is larger in the pure 
p-mode region. This causes the departure from a linear splitting 
due to cubic order effects including near degeneracy to be of 
the same order of magnitude as the latitudinal shear contribution 
(Fig.©. 
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Fig. 7. Sf, m (pi) versus Sf, m (g2) due to cubic order effects and 
near degeneracy (purple), and to a latitudinal shear (green), 
for ( — 1 modes in an evolved 8.5 M model (see Table |2). 
The rotation rate increases from left to right: from 18% fit to 
34%Qfr. The er r or bar s for V 1449 Aquilae have been taken from 
iDegroote et"al] (120091) . 



6.2. Departure from a linear splitting as a function of the 
rotation rate 

We studied the relative values of third order, near degeneracy 
and latitudinal shear contributions to the splitting as a function 
of the rotation rate mode by mode. Figure|6]shows three different 
cases: 

- the g5 mode illustrates the the case of high order g-modes (or 
pure g-modes) where the cubic order contribution overtakes 
the latitudinal shear contribution as soon as the mean rota- 
tion rate exceeds 15 % Q* (which corresponds to a rotation 
velocity of around 80 km.s 1 ). 

- the g2 mode is located in the mixed mode region where cubic 
order effects are larger than latitudinal differential rotation 
ones. 

- the pi mode is still a mixed mode but with a nature closer 
to a pure p-mode for which near degeneracy is no longer 
negligible. As a result, the total cubic contribution including 
degeneracy is of same order as a latitudinal shear of 12.7 %. 

6.3. How do we disentangle the two effects? 

In the previous subsection, we have shown that for a massive star 
on the main sequence, in the frequency range where we expect 
pulsation modes, it is not easy to conclude whether a departure 
from a linear splitting is due to latitudinally differential rotation 
or to cubic order effects. Here we suggest a method to disentan- 
gle the two effects. 

Let us consider two modes: g2 and p\ as representative of 
a mixed mode close to g-modes and a mixed mode close to p- 
modes, respectively, as seen in the previous section. In Fig. [7j 
S i (p i ) versus Si(g2) is plotted for the two different assumptions 
(latitudinal shear or cubic order effects with near degeneracy cor- 
rections) for different rotation rates (Q < 35%Q^). Note that 
for g-modes, f3 is roughly equal to 1 /2, and for p-modes /3 ap- 
proaches 1. Accordingly, the ratio between S' at (pi) and Sj at (g2) 
does not depend on AQ and is approximatively equal to 2. On 
the other hand, the curve S i (p i ) as a function of S i (g2) for split- 
tings including cubic order effects with near degeneracy correc- 



tions starts to deviate from a slope of 2 when the rotation rate is 
large enough. The deviation increases with the rotation rate as 
expected. 

Let us now define S° bs as the observed splittings for / = 1 
modes for a fast rotator with a uniform rotation profile Q[["'\ 

Their values are then assumed to be given by S° bs = S° (^^ ue ) 
(Eq.(l22li) i.e. rotational splittings computed up to cubic order, 
including near degeneracy corrections. 

If one misinterprets S° bs as due to a latitudinal shear, S° bs is 
assumed to obey Eg. (132b. The point representing S* s (Pi) as a 
function of S° bs (g2) ought to lie on the straigth line with slope 2 
in Fig.|7j We represent this point with the splittings S* s (pi) and 
S* s (g2) computed according to Eq.(l22l assuming a rotation rate 
of 15%Qfc, v = SO km.s' 1 . As the observed point is not located 
on the straight line given by S^'(pi) = 2S[ a 1 '(g2), we are then able 
to conclude that the deviation from a linear splitting is not due 
to a latitudinal shear. 

In the case of massive stars on the main sequence, the ef- 
fects of latitudinal differential rotation and cubic order or near 
degeneracy are of the same order of magnitude in the frequency 
domain where we expect to observe oscillation modes. This is 
mostly due to the mixed nature of modes around the fundamen- 
tal frequency, and should therefore depend on the evolutionary 
stage of the star. But if two individual rotational splittings are 
available (one g-mode under the mixed mode region, and one 
p-mode above it), this method helps to distentangle whether a 
departure from a linear splitting is due to cubic order contribu- 
tions including accidental degeneracy, or latitudinal shear. 



7. The case of solar type stars 

Low mass main sequence stars are known to be slow ro- 
tators. Indeed due to their outer convection zone, they un - 
dergo magnetic braking during their evolution dKawalerl 119881) . 
Observational evidence exists for surface latitudinal differential 
rotation (coming from stellar spots due again to their outer con- 
vection zone). Hence, for these stars, the averaged surface rota- 
tion rate Q is much smaller compared to that of more massive 
stars, such as /3 Cephei stars, and AQ = Q. equa - Q. po ie, the dif- 
ference between the rotation rates at the equator and the poles, 
can be large ( 25% for the Sun, betw een 1% and 45% for a star 
like Procyon, iBonanno et all 120071) . One therefore expects that 
latitudinal corrections to the splittings will dominate over cubic 
order effects which are negligible. As illustrative examples, we 
studied the case of HD 181906 and HD 181420, which are two 
solar like stars observed by CoRoT during t he first long r u n, and 
whi ch lightcurves have been analysed by iGarcia et al.l Q2009) 
and lBarban et al.l d2009i) respectively. 

Before comparing the different contributions to the rotational 
splitting, one should wonder whether a perturbative approach up 
to the cubic order accounting for near degeneracy effects is valid 
for seismic interpretation purpose for this type of st ar. To answer 
this is sue, we rely on the validity study done in ISuarez et al.l 
(l2010l) . This study has been performed for a poly tropic model 
of 1.3 M , that can be representative of HD 181420. Therefore 
we consider this study as well adapted in order to determine the 
validity of our approach for computing high order pressure pul- 
sation modes which are propagating in the o uter layers of a solar 
like star. In this work (ISuarez et al.l 1201 Ol) show that for rota- 
tional velocities under approximately 40km. s , the structures 
of the frequency spectra computed with a non perturbative and 
their perturbative method are very similar. Here, we study the ro- 
tational splitting (that is a frequency difference which removes 
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part of second order effects), with a cubic order perturbative ap- 
proach accounting for the effect of near degeneracy. We then 
consider our approach as valid for rotation velocities at stake in 
the stars we study in this article. 



7. 1 . Competition between the three effects 

Here we investigate the order of magnitude of deviations from 
Eq. (f24b for an M — 1.36 M , R — 1.63 R main sequence stellar 
model (see Table [3] for the stellar parameters taken as a model 
representative of HD 181420). 



Table 3. Stellar parameters of the model of HD 18 1420 (Sect.6) 



M = 1.36M 




R = 1.63 Rq 


L = 4.4 L Q 




X = 0.7 


P c = 183dyn.cirT 2 


Pc 


= 5.1xl(T 10 g.cnr 3 



In Fig. [8] we plot the three kinds of contributions - i.e. cubic 
order, cubic order with near degeneracy, and latitudinal shear - 
to (€ — 1, ±1) splittings as a function of the central mode fre- 
quency for two different rotational angular velocities (3.7% Q* 
and 7.4% fi^, which correspond to 15km. s _1 and 30km. s -1 , re- 
spectively, for the model described in Table©. Once again, these 
contributions show a peculiar behavior in the mixed mode re- 
gion, where the two types of contributions can be of the same 
order. We focus here on higher frequencies, since the oscilla- 
tion modes of HP 181420 are det ected in the 1.5 - 2 mHz fre- 
quency range (iBarban et al.L 12009). Therefore, these plots show 
that even for rotation rates which are high for this type of star 
(7.4% fi*), cubic order contributions (with or without near de- 
generacy) can be neglected compared to additional terms due to 
latitudinal shear. This leads to the conclusion that for HD 1 8 1420 
in particular - as we expect for solar-like stars in general - the 
cubic order contributions to the splitting can safely be neglected 
in comparison with potential latitudinal shear contributions. The 
same computations have been performed for a 1.2M 1.4 Rq 
stellar model representative of HD 181906, and lead to the same 
conclusions. 



7.2. Determination of a latitudinal shear from seismic 
observations: 



The re sults of the seismic an alysis provided in IBarban et al.l 
d2009t) and lGarcfa et all (120091) both find significantly different 
values for the low frequency peak in the Fourier spectrum and 
what is interpreted as the rotational splitting. If a uniform sur- 
face rotation is assumed, it is not possible to reproduce such dif- 
ferences between these observables. 
Let us assume a rotation profile of the form: 



fi(<9) = fi o -Aficos 2 



(39) 



Where fin is the rotation surface angular velocity at the equator, 



fie, 



quotor: 



and Afi = fi e , 



then given by : 



quotor 



Qpoie- The rotation frequency is 



v rot (0) 



"eq 



Afi 

fin" 



COS" 



(40) 



where v eq and v ro , correspond to the equatorial rotation rate (i.e. 
fin) and the rotation frequency at the colatitude 9 in //Hz respec- 
tively. 

Using again the rotation law Eq.(l39l to derive the rotational 
splitting from Eq.(f27b. we obtain for ( = l,m = ±1: 



lat fin / 1 Afi 
Sl - 1- fi^^Sfin" 



which can be written under the following form: 

Afi 

fiu" 



Vsplit = V eq j6 1 



(41) 



(42) 



where v sp iu corresponds to Sij. 

Then Eq.(l40b and Eq.(l42b are two equations with three un- 
knowns: v eq , cos9 and Afi/fin. Dividing Eq. d42l by Eq.d40l) 
leads to : 



Afi 
fio" 



1 -/3d 



cos#- 



,Bd 



Where we have introduced the parameter: 

, Vrot 

d - 

Vsplit 



(43) 



(44) 
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q 

< 




(pd/5) 



cos 9 



Fig. 9. Possible values for AQ/Clo, depending on the parameter 
d Eq.(l44l and the value of the integral /3 Eq. (IB.5l ). The line col- 
ors corresponds to the value of (d x /?) which ranges here from 
0.3 (in blue) to 1.5 (in red). In grey is the acceptable domain of 
values for the latitudinal gradient (Eq. [46}. This plot has been 
performed for high order p modes. 



We then determine AQ/Qo as a function of cos 9, along with the 
following constraints: 



Considering the uncertainties attached to the seismic observ- 
ables, the result obtained for the latitudinal shear is presented in 
Tab. [5] These results show that the hypothesis of latitudinal dif- 
ferential rotation is consistent with the available observables . 
Therefore, the latitudinal differential rotation profile is a reliable 
assumption concerning HD 181906. Moreover, The results de- 
fine a range of possible values for different characteristics of the 
rotation profile: the latitudinal shear, the rotation surface velocity 
at the equator, as well the inclination angle. The correspondant 
error bars are discussed in Sect. 17.51 Both scenarii give the same 
conclusionsvery similar results. 

Table 5. Latitudinal shear, equatorial velocity and inclination an- 
gle obtained for HD 181 906, using models with rotation uniform 
in depth. The two last lines correspond to i - arcsin(vsin!'/v e? ) 

with v sin/ = ( 10+ l)km.s~' (iBrunttl , l2009h . or with vsiru = 
(16+1) km.s 1 dNordstrom et all 120041) . Note that the two spots 
give hardly the same results, we present here mean values com- 
puted for the two surface rotation velocities. 





Scenario A 


Scenario B 


An/fio = 

v er? (km.s- 1 ) = 

i (Brunt 2009) (°) = 

i (Nordstrom 2004) (°) = 


0.55 ± 0.21 

41.3 ± 2.5 

26 ± 2 

23 ± 3 


0.59 ± 0.18 

43.8 ± 1.6 

24 ± 2 

21 ± 2 



n 
° ; 2 



cos0 € [0;1] 

AD, 

6 [-1;1] 

Oo 



(45) 
(46) 



As illustrated in Fig. [9] the latitudinal shear is a hyperbolic 
function of cos 9, centered in J/M/5. Depending of the values 
of the product /3d of the observationnal parameter d and of the 
integral /?, graphically we find different possible ranges of values 
for Af2/f2o and cos 9. 

7.3. The case of HD 181906 

HD 181906 is an F8 dwa rf for which f undamental parameters 
have been determined by iBrunta (2009), considering the pres- 
ence of a background star : m v - 7.65, L/L Q = 3.32 + 0.45, 
T eff = (6300+150)K,[Fe/H] = -0.11 ±0.14, M/M Q = 1.144 + 
0.119, vsini = (10± l)km.s _1 . Note that iNordstrom et all 
(2004) had found vsini = (16±l)km.s~ 1 considering HD 
181906 as single, with no background star. It has been observed 
during the fi rst long run of CoR oT, and its light curve has been 
analysed by IGarcia et al.l d2009l) . They find two possible inter- 
pretations for the mode identification, as listed in Table |4] 



Table 4. Results con cerning rotation from the analysis of HD 
181906 performed in IGarcia et al.l (120091) . v lot i and v rot 2 stand 
for the two low frequency peaks attributed to the rotation of two 
different spots on the star surface. These values correspond to a 
product (3d ranging from 0.65 + 0.03 to 0.70 ± 0.03. 



Scenario A 



Scenario B 



Vrotl 0"Hz) 
v rot2 (//Hz) 
v sp ii, QuHz) 



4.0 : 

4.45 



5.8 ± 0.14 



0.15 

: 0.15 

6.1 



7.4. Solar or anti-solar latitudinal shear for HD 181420? 

HD 181420 is an F2 main sequenc e star whose fu ndamental pa- 
rameters have been determined by |Bruntt (2009) : m v = 6.57, 
L/L = 4.28+0.28, T eff = (6580±105)K, [Fe/H] = 0.00+0.06, 
M/M = 1.31 ±0.063, vsini = (18 ± ^km.s 1 . Ithas been ob- 
served during the first long run of CoR oT, and its light curve has 
been analysed by iBarban et al.l (12009b . They find two possible 
interpretations for the mode identification, as listed in Table [6] 
Later Bedding & Kjeldsen (2010) used some empirical scaling 
method that seemed to favor scenario 1 . 



Table 6. Results co ncerning rotation fro m the analysis of HD 
1 8 1420 performed in IBarban et al.l (120091) . 



Scenario 1 



Scenario 2 



v rot (^Hz) : 
Vspiii (y"Hz) 



4.5 



2.59 ± 0.38 



1.5 
3.29 



0.17 



Concerning HD 1 8 1420, our model gives unexpected results: 
in both scenarii, the latitudinal shear is found negative, i.e. the 
pole rotate faster than the equator, which is the opposite of what 
is known for the Sun. In order to apprec iate the reliability o f 
this result, one should refer to the work of iKapvla et al.l (1201 lb . 
who study the impact of rotation on turbulent angular momen- 
tum and heat transport in solar like stars convective zone, by the 
mean of direct numerical simulations of turbulent convection in 
spherical geometry. According to the authors, the rotation pro- 
file varies from anti-solar (equator rotates slower than poles) for 
low Coriolis number to solar (equator rotates faster), with a tran- 
sition around Co = 3. The Coriolis number which measures the 
impact of rotation on turbulent motion is given by: 



0.14 



Co = 2 fio T c 



(47) 
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Where r c is the convective overturn ing time scale -computed in 
the stellar evolution code CESAM dMorel & Lebretoni l2008h -. 
and Qo is taken as the angular velocity at the equator. For HD 
181420, when rotation is supposed to be uniform in depth, the 
Coriolis number is (3.2 + 0.5) for scenario 1 and (4.0 + 0.3) 
, i.e. both values are only slightly above the transition thresh- 
old between anti-so lar and solar surface rotation (see Fig. 17 of 
Kap yla et all 1201 it) . Therefore, it is difficult to confirm that the 
physics at stake in the convective envelope of HD 181420 can 
lead to anti-solar rotation profile for HD 181420. 

Let us refine the modeling of HD 181420, taking a 
tachocline-like profile. As explained in Sect. 15.21 it consists in 
assuming uniform rotation (Q = Qq) in the inner layers, and 
latitudinal differential rotation (Q(0)) in the convective envelop. 
The rotational splitting given by Eq.fiTTi is no longer relevant, 
and should be replaced by Eq. (l3"7T i. 

In this case, we found different values for fio than in the 
former case, as a consequence, the Coriolis number reaches 
(8.9 + 1.7) for scenario 1 and (11.8 + 0.7) for scenario 2, i.e. 
high above the transition limit between anti-solar and solar ro- 
tational shear. The integral a reaches a constant value of 0.45 
for high order p-modes, which changes the relationship between 
the colatitude of the spot and the latitudinal shear ( 1431 . and as a 
consequence changes the domain of possible shear. Positive lat- 
itudinal shear corresponding to solar type latitudinal rotation are 
obtained for the two scenarii, which is this time fully consistent 
with the Coriolis number values. The results obtained with this 
tachocline-like model are listed in Tab. [7] 



Table 7. Latitudinal shear, equatorial velocity and inclina- 
tion angle obtained for HD 181420, using models with 
tachocline-like rotation profile. The last line corresponds to i - 
arcsin (v sin i/v eq ) with v sin; = (18+ l)km.s _1 dBrunttl l2009h 





Scenario 1 


Scenario 2 


Afl/n = 

v eq (jMz) = 

\ eq (km.s- 1 ) = 

i (Brunt 2009) (°) = 


0.50 ± 0.45 

1.0 ± 0.2 

45.6 ± 8.7 

23.2 ± 7.4 


0.66 ±0.11 
1.3 ± 0.1 
60.1 ± 3.5 
17.4 ± 2.0 



The case of HD 181420 is particularly interesting as a sim- 
ple model of rotation uniform in depth and differential in lati- 
tude leads to non-physical latitudinal shear. Only a little more 
sophisticated model, where rotation varies in depth according 
to a tachocline-like profile, is fully coherent with the physics 
of the convective zone as well as the observables. We also give 
a range of possible values for the rotation profile parameters: 
the latitudinal shear in the convective zone, the rotation surface 
velocity at the equator and the incination angle. To conclude, 
for this star, we find that rotation profile inside the star should 
rather be a tachocline like profile, with a solar type latitudinal 
rotation. Moreover, these results, when compared to those of 
iMosser et al.l d2009l) . seem to favor scenario 1 since the rota- 
tional frequency at the equat or in the hypothesis of scenario 1 
is closest to the one found in IMosser et al.l d2009l) by spot mod- 
elling - (5.14 ± 0.07) x lO^rad.s- 1 -. 

7.5. Discussion 

As already menti onned in the beginning of this section, 
ISuarez et al.l (1201 Ol) have found that second order perturbative 
methods including near-degenaracy corrections are valid for ro- 



tation velocity up to approximately 40km. s . In this study, third 
order perturbative methods have been used for the computation 
of rotational splittings in Sect. 7.1. In Sect. 7.2, is given a sim- 
ple recipe which allows to compute the latitudinal shear given 
seismic observables. This formulation rely on the validity of the 
third order perturbative method, but the only quantity provided 
by the modeling are/? (Eq. IB.5l ) and/3 cz , i.e. quantities computed 
from first order eigenfunctions. The question is then: are eigen- 
functions suffitiently varing with rotation rate (under 60km. s" 1 
i.e. 15 % Qj), to impact the value of/?? After verifications, it ap- 
pears that for high order acoustic modes in moderately rotating 
stars (under 15%Q/t)/3 - 1, and/3 cz ^ 0.45 are satisfying evalu- 
ations. Therefore the recipe is valid for rotation rates at stake in 
Sect. 7.3 and 7.4. Note the convenience of the proposed recipe 
which only need for seismic observables, and the values of (3. 

The large error bars obtained in particular for the values of 
latitudinal shears (Tab. [5][7j is not only due to the observational 
uncertainties, but it can also be attributed to the simplicity of the 
spot model used. This model does not account for spot parame- 
ters such as for example a spot lifetime, or the spot distribution 
on the observed stellar disc (for more sophisticated modelling 
see IMosser et all 120091) . We only consider a unique spot of in- 
finite lifetime. Note that the use of a more complicated model 
would require more observational constraints than only one spot 
rotation signature. When only mean values of rotationnal split- 
tings are available, and the observational error ba rs on the low 
frequ ency signature of a spot rotation are large dBarban et all 
120091) . we are able to give general conclusions concerning the 
rotation profile -i.e. uniform in depth or tachocline like, solar or 
anti-solar-, but no reliable numerical values caracterising it. 



8. Conclusions 

With the help o f the perturbative approach established in 
ISoufi et akl dl998l) . we investigated second and third order con- 
tributions of the Coriolis and the centrifugal accelerations to p 
and g-modes frequencies, as well as near degeneracy effects on 
the rotational splittings. Their effects were then compared with 
those of a latitudinal shear. We studied two types of oscillating 
stars. 

For an evolved model of a j3 Cephei, the effects of near de- 
generacy, cubic order perturbations and latitudinal shear are of 
the same order of magnitude in the frequency range, relevant to 
such stars - i.e. low order p- and g-modes - and for reasonable 
rotation rates ranging from 15% Q* to 30%^. Nevertheless, 
when two individual splittings for modes of a different nature (a 
pure g and a mixed mode for example) are available, a method is 
proposed to distinguish between a latitudinal shear and the other 
effects. 

For solar-like stars such as HD 181420 and HD 181906, 
which are mostly moderate rotators and oscillate with high order 
p modes, cubic order effects on frequency splittings are shown 
to be small in front of the effects of latitudinally differential rota- 
tion. Therefore, given a splitting and a rotation period signature, 
it is possible to infer a range for the latitudinal shear coefficient 
AO/Qo- Although no precise determination of latitudinal shear 
is possible unless the spot latitude is fully determined, we have 
been able to determine the most reliable rotation profile for each 
of the two stars. More over the determination of the latitudinal 
shear by our seismic method can be taken as a constraint to be 
added to other observational constra ints, such as those coming 
from spectropo larimetric results (e.g . iDonati et all 1201 Oh or spot 
modelling (e.g. IMosser et all 120091) . 
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Appendix A: Near degeneracy corrections to the 
splittings of high order p-modes 

This section is dedicated to a qualitative estimation of the near 
degeneracy corrections to rotational splittings: 



;.deg 
5 f=l,m=l 



^ 



+ <H} 



V A -i + tt-i 



J f=l,m=l 



- -(A, - A_0 



(A.l) 
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where Ai and Hi are normalized by the break-up frequency Q^. 
The oscillation frequency given in Eq.© can be rewritten in the 
following the generic form : 

o"„/, m = <r 0n{ + m—f3„j + m dn,f,\ m \ + mt n j,\ m \ 

where d, h t,\ m \ and f B /,|m| represent the second and third order con- 
tributions. Therefore the term Ai and A_i appearing in Eq. (l22t . 
for the coupling of an I — \,n and an I — 3,«' mode, can be 
written as: 

n & 

Ai = A u + —(/? n ,i -y3„', 3 ) + (rf«,i,i - rf»',3,i) + (f«,i,i - f«',3,i) 

(A.2) 
where A stands for cr ( l , — cr5 , ,. the difference between 

U,»,l 0,n ,j 

the eigenfrequencies without rotation. 

For high frequency p-modes, the radial contributions are 
concentrated in the outer layers and are then nearly the same 
for («, I — 1) and («', I - 3). It is then legitimate to neglect the 
differences in their radial contributions. We then omit the radial 
order subscripts n, n' . For the same reason, the Ledoux constants 
are also quite similar/?! ~ /?3- Then 

A! * A + (d u - d xx ) + U - f 3 ,i) 

Similarly 

A_i * A + {d„ X \ -dxi) - (fi,i -*3,i) 

Hence, the quantity (Ai - A_i) in Eq. dA.ll ): 

A_i - A! ~2(f u -f 3 ,i) (A.3) 

is of third order. One can then approximate A_i ~ Ai in the 
square root of Eq. dA.U . and the correction due to near degener- 
acy is given by: 

Q deg Q ND ^ , t t \ 

i> e= i - jg = i ~ (n,i -*3,i.) 



V A ? + ^ 



a/ a ? + ^ 2 i 



(A.4) 

For the coupling term, one can separate the second Hi and third 
'Hi order terms as H m = <H% + mH?,. Thus 'Hi = H 2 + < H 3 , 
and H-i - H 2 - "K 3 . With |*H 3 | « \<H%\, the near degeneracy 
correction to the splitting becomes: 



jdeg 
'f=l 



>C=\ 



(tl,l-t 3 ,i) + 






(A.5) 



which shows that for slow rotators, such as HD 181420, if cubic 
order effects are neglected, then the near degeneracy contribu- 
tion is zero, and the splitting is a linear function of rotation. If 
cubic order effects are included, then near degeneracy correc- 
tions affect the splitting, and the departure from a linear splitting 
is dominated by distorsion (predominantly in "Tr^) . 
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a. for£= l,m = ±1 (i.e. A = 2): 



Appendix B: Contribution of the latitdinal shear to 
the splittings 

In order to be able to compute the splittings from Eq.([T} and 
Eq.dT7t. one must specify a rotation law. It is convenient to as- 
sume the following form: 



M.l 



Qo 

-Qk 



P 



1 



1 AQ 

5 "o7 



b. for ^ = 2, m= +1 (i.e. A = 6) : 



Q(r, G) = V Q 2s (r) (cos 0) 2s 



(B.l) 



s=0 



S2J = * 



1 AQ 
/3 + - — (8 7-3,6) 
/ 12 



where 6 is the colatitude and we have limited our investigation 
Smax - 1 ■ The surface rotation rate at the equator is Q(r = R, — 
nil) = n {T = R). 

Inserting Eq.(l25Tl into Eq.(IT7b yields the following expres- 
sion for the generalized splitting: 



c. for t = 2, m - +2 (i.e. A = 6) : 



S2 ' 2 = a 



1 AQ 
2/3 -- — {2y+ P) 

I kl 



(B.13) 



(B.14) 



(B.15) 



-R 



1 



s=2 



Q* Jo "k ^ 



Appendix C: Conditions for significant near 
degeneracy 



with (see lGoupilil201lh . 

H m , s (Q) = -1/1 [R. (^ 2 - 2£& + £ 2 (A - 1)) + Q s ^)]p r 2 dr 

where R s and Q s depend on Q 2 , Q4 and A = ({€ +1). 

For example, if shellular rotation is expected, then Q(r, 9) - 



As already mentionned, near degeneracy between two modes oc- 
curs whenever their frequencies are close, their azimuthal num- 
bers are equal, and their angular degrees differ by 2. However the 
magnitude of the near-degenerate corrections to both frequen- 
cies also depends on the magnitude of the coupling term H m . In 
turn the magnitude of H m depends on the nature of the involved 
modes whether they are g-modes, p-modes or mixed modes with 
a dominant p or g nature. 

The g-mode spectrum is much denser than the p-mode one. 
Hence, as shown in Fig. C.2, A,„ is much smaller than for p- 
modes. This ought to favor near degeneracy. However, the cou- 
pling term H„, for g-modes is much smaller than for p-modes. 
If we consider latitudinally differential rotation only, Q 2 ;,j = This is due tQ ^ fact ^ distorsion effects ^ small for g . 

0, 2 are depth independent, Qo = Q eq uator, ^2 = -AQ, and Q4 = modes and therefore the overall (second and third order) con- 



Qo(r) and s max = 0. Moreover, Q2, Q4 
forj = 0,2, and: 



0, i.e. R 



Q/. 



[g-2^ h + Af h ]p r 2 dr 



Qj = 
(B.3) 



0. R s and Q s are constant and: 



Q 1 s=1 

k k I^o 



tribution to H m remains small. As a result, the coupling term is 
much smaller than A,„ for g-modes, which are then hardly cou- 
pled. 



with/5 and y defined as: 



P = \ J [£- 2£ft + (A - D4 2 ] Por 2 dr 



y 



"I Jo 



S PoT dr 



and I being the inertia of the mode: 

1= f (^ 2 +A4 2 )p r 2 dr 
Jo 

Then the rotational splitting of a {n, 6, m} mode is given by 
S e>m = -j- [Q 0y 8 + R yS + Q r + m 2 (R^ + Qf r )] (B.8) 



(B.5) 
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In lGoupill d201lh . Q s and R s are given for s = 0, 1 , 2: 

. AQ 2A - 1 

Rn 



R 



Q 



Qo 


4A- 


3 


AQ 


2 




Q 4A-3 




AQ 


2(1- 


3A) 


Q 


4A 


-3 


AQ 


10 





1 Qo 4A-3 

For £ = 1 and t — 2 modes, this yields: 



Fig. C.l. Inertia of axisymmetric I = 1 modes as a function of 

to the radial order, n, for an evolved 8.5 M 0y 8 Cephei model, uni- 

(B.9) formly rotating at 30% Q A ,i.e. 160/tm.s^ 1 (see Tabled Sect.4.1). 

(B.10) 

(B.ll) 

(B.12) 
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Fig. C.2. The terms A m and 1ri m which are involved in near de- 
generacy corrections (see Eq.tfTJt) for (. — 1 modes computed 
for the same f3 Cephei model as in Fig.C.l. A logscale is used 
for the y axis. 
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Appendix D: Values for third order contributions to the frequencies and to the splittings 



Table D.l. Scaled contributions to the splittings of t — 1 g-modes due to: implicit third order terms in the eigenfrequency (T e,gm ), 
Coriolis effects (T T ), distorsion (T D ), and coupling of the two (T c ). The impact of near degeneracy is so small that its contribution is 
fully negligible and is therefore not listed. The first row lists the radial order and the second row the centroid mode m = frequency. 



n 


a-o 


rpeigen 
1 3 


n 


1 3 


7* 


-8 


0.25 


0.612E + 02 


0.529£ + 04 


-0.390E + 01 


-0.557 E + 03 


-7 


0.29 


0.705£ + 02 


0.391£ + 04 


0.214£ + 02 


-0A24E + 03 


-6 


0.34 


0.770£ + 02 


0.285£ + 04 


0.379£ + 02 


-0.362E + 03 


-5 


0.41 


0.943£ + 02 


0.208£ + 04 


0.518£ + 02 


-0.354E + 03 


-4 


0.52 


0.123£ + 03 


0.153£ + 04 


0.683£ + 02 


-0.371£ + 03 


-3 


0.70 


0.155£ + 03 


0.111E + 04 


0.905£ + 02 


-0.374£ + 03 


-2 


1.06 


0.850£ + 02 


0.776£ + 03 


0.116£ + 03 


-0.318£ + 03 


-1 


3.52 


-0.269£ + 04 


0.833£ + 03 


-0.439£ + 03 


0.851£ + 03 



Table D.2. Scaled contributions to the splittings of ( = 1 p-modes due to third order effects divided by the square of the central 
mode frequency. The first row lists the radial order and the second row the centroid mode m - frequency. 



n 


0"o 


rr'K 


Tlloi 


T?/o* 


m< 


T?K 


-1 


3.52 


-0.211 E + 03 


0.612E + 02 


-0355E + 02 


0.6S1E + 02 


-0.782£ - 05 


1 


4.73 


-0.140E + 03 


0.409£ + 02 


-0.181£ + 02 


0.433£ + 02 


0.691 E + 00 


2 


5.78 


-0.106£ + 03 


0.285£ + 02 


-0.115£ + 02 


0.343£ + 02 


0.174£ + 02 


3 


6.85 


-0.862£ + 02 


0.181£ + 02 


-0.554£ + 01 


0.264£ + 02 


0.318£ + 02 


4 


7.96 


-0.121 E + 02 


0.112£ + 02 


-0.666£ + 00 


0.192E + 02 


0.380£ + 02 


5 


9.11 


-0.636£ + 02 


0.767£ + 01 


0.217£ + 01 


0.145£ + 02 


0.435£ + 02 


6 


10.30 


-0.566£ + 02 


0.581£ + 01 


0.384£ + 01 


O.llOf + 02 


0.447£ + 02 


7 


11.49 


-0.514£ + 02 


0.494£ + 01 


0.456£ + 01 


0.884E + 01 


0.472£ + 02 


8 


12.67 


-0.477£ + 02 


0.440£ + 01 


0.472£ + 01 


0.762E + 01 


0.503£ + 02 


9 


13.82 


-0.447£ + 02 


0.374£ + 01 


0.471£ + 01 


0.688E + 01 


0.540£ + 02 


10 


14.96 


-0.417E + 02 


0.290£ + 01 


0.489£ + 01 


0.601£ + 01 


0.555£ + 02 


11 


16.11 


-0.387£ + 02 


0.221£ + 01 


0.523£ + 01 


0.494E + 01 


0.543£ + 02 


12 


17.28 


-0.360£ + 02 


0.181£ + 01 


0.552£ + 01 


0.395E + 01 


0.521 E + 02 


13 


18.46 


-0.338£ + 02 


0.1615 + 01 


0.570£ + 01 


0.316E + 01 


0.516E + 02 


14 


19.63 


-0.318£ + 02 


0.149E + 01 


0.580£ + 01 


0.252E + 01 


0.50SE + 02 


15 


20.81 


-0.301£ + 02 


0.140£ + 01 


0.585£ + 01 


0.196£ + 01 


0.498£ + 02 



16 



R-M. Ouazzani, and M-J. Goupil: Rotational splittings for slow to moderate rotators: 
Table D.3. Different order contributions to the mode frequencies for various radial orders, n. All contributions are scaled by Q^. 



n 


<ro 


<T\ 


eigen 

u n 




a? 




*? 




-8 


0.25 


0.777£-01 


0.9895 - 


-02 


0.268£ - 


-01 


-0.1675- 


-03 


-7 


0.29 


0.782E-01 


0.S57E - 


-02 


0.2385 - 


-01 


-0.1895- 


-03 


-6 


0.34 


0.7905-01 


0.7135 - 


-02 


0.2065 - 


-01 


-0.2285 - 


-03 


-5 


0.41 


0.802£-01 


0.562£ - 


-02 


0.1735- 


-01 


-0.2935 - 


-03 


-4 


0.52 


0.8195-01 


0.4075 - 


-02 


0.1415- 


-01 


-0.3925 - 


-03 


-3 


0.70 


0.8395-01 


0.252£ - 


-02 


0.1085- 


-01 


-0.5295 - 


-03 


-2 


1.06 


0.8595-01 


0.954£ - 


-03 


0.7375 - 


-02 


-0.7295 - 


-03 


-1 


3.52 


0.142£ + 00 


0.677£ - 


-02 


0.3255 - 


-02 


-0.1015- 


-01 


1 


4.73 


0.141£ + 00 


0.497£ - 


-02 


0.2475 - 


-02 


-0.1435- 


-01 


2 


5.78 


0.1415 + 00 


0.4015 - 


-02 


0.2055 - 


-02 


-0.1885- 


-01 


3 


6.85 


0.142£ + 00 


0.3375 - 


-02 


0.1765- 


-02 


-0.2245 - 


-01 


4 


7.96 


0.142£ + 00 


0.2885 - 


-02 


0.1535- 


-02 


-0.2535 - 


-01 


5 


9.11 


0.1435 + 00 


0.25 15- 


-02 


0.1355- 


-02 


-0.2865 - 


-01 


6 


10.30 


0.144£ + 00 


0.2215 - 


-02 


0.1215- 


-02 


-0.3195- 


-01 


7 


11.49 


0.1455 + 00 


0.1975- 


-02 


0.1095- 


-02 


-0.3565 - 


-01 


8 


12.67 


0.1455 + 00 


0.1785- 


-02 


0.9975 - 


-03 


-0.3975 - 


-01 


9 


13.82 


0.1465 + 00 


0.1635- 


-02 


0.9195- 


-03 


-0.4385 - 


-01 


10 


14.96 


0.1465 + 00 


0.1505- 


-02 


0.8525 - 


-03 


-0.4715- 


-01 


11 


16.11 


0.1465 + 00 


0.1395- 


-02 


0.7945 - 


-03 


-0.4995 - 


-01 


12 


17.28 


0.147£ + 00 


0.1295- 


-02 


0.7435 - 


-03 


-0.5265 - 


-01 


13 


18.46 


0.147£ + 00 


0.1215- 


-02 


0.6975 - 


-03 


-0.5575 - 


-01 


14 


19.63 


0.1475 + 00 


0.1145- 


-02 


0.6575 - 


-03 


-0.5885- 


-01 


15 


20.81 


0.147E + 00 


0.1075- 


-02 


0.6215- 


-03 


-0.6195- 


-01 



n 


of m 




T 

^3 




of 




of 




o- dcg - o-' 


lodeg 


-8 


0.1865- 


■04 


0.1935- 


-02 


0.1995- 


-04 


-0.5285 


-03 


< 10" 


5 


-7 


0.2595 - 


04 


0.1725- 


-02 


0.3255 - 


-04 


-0.4395 


-03 


< 10" 


5 


-6 


0.1855- 


04 


0.1515- 


-02 


0.4235 - 


-04 


-0.3945 


-03 


< 10" 


5 


-5 


0.5825 - 


05 


0.1305- 


-02 


0.5155- 


-04 


-0.3785 


-03 


< 10" 


5 


-4 


-0.2055 ■ 


-04 


0.1055- 


-02 


0.6075 - 


-04 


-0.3625 


-03 


< 10" 


5 


-3 


-0.6965 


-04 


0.7675 - 


-03 


0.6775 - 


-04 


-0.3135 


-03 


< 10" 


5 


-2 


-0.1785 


-03 


0.4295 - 


-03 


0.6405 - 


-04 


-0.2075 


-03 


< 10" 


5 


-1 


0.1395- 


03 


0.1735- 


-03 


-0.9605 


-04 


0.1875- 


-03 


0.5075 - 


-03 


1 


0.8835 - 


04 


0.1445- 


-03 


-0.6555 


-04 


0.1565- 


-03 


0.1005- 


-02 


2 


0.6355 - 


04 


0.1235- 


-03 


-0.5065 


-04 


0.1515- 


-03 


0.2625 - 


-02 


3 


0.3525 - 


04 


0.9305 - 


-04 


-0.2885 


-04 


0.1375- 


■03 


0.4475 - 


-02 


4 


0.4965 - 


05 


0.6705 - 


-04 


-0.4025 


-05 


0.1165- 


-03 


0.5965 - 


-02 


5 


-0.1925 


-04 


0.5265 - 


-04 


0.1495- 


-04 


0.9985 - 


-04 


0.2035 - 


-02 


6 


-0.3595 


-04 


0.4505 - 


-04 


0.2995 - 


-04 


0.8545 - 


-04 


0.2995 - 


-02 


7 


-0.4895 


-04 


0.4285 - 


-04 


0.3955 - 


-04 


0.7675 - 


-04 


0.4155- 


-02 


8 


-0.5955 ■ 


-04 


0.4205 - 


-04 


0.4525 - 


-04 


0.7295 - 


-04 


0.1875- 


-01 


9 


-0.7135 


-04 


0.3905 - 


-04 


0.4925 - 


-04 


0.7185- 


-04 


0.2445 - 


-01 


10 


-0.8305 


-04 


0.3275 - 


-04 


0.5535 - 


-04 


0.6795 - 


-04 


0.2955 - 


-01 


11 


-0.9065 


-04 


0.2695 - 


-04 


0.6375 - 


-04 


0.6015- 


-04 


0.3365 - 


-01 


12 


-0.9645 


-04 


0.2375 - 


-04 


0.7215- 


-04 


0.5165- 


-04 


0.3795 - 


-01 


13 


-0.1015 


-03 


0.2255 - 


-04 


0.7965 - 


-04 


0.4415 - 


-04 


0.4345 - 


-01 


14 


-0.1055 


-03 


0.2225 - 


-04 


0.8625 - 


-04 


0.3745 - 


-04 


0.4975 - 


-01 


15 


-0.1105 


-03 


0.2215- 


-04 


0.9225 - 


-04 


0.3095 - 


-04 


0.5645 - 


-01 
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